Off-diagonal Long-Range Order and Supersolidity in a Quantum Solid with Vacancies 
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We consider a lattice of bosonic atoms, whose number A'^ may be smaller than the number of 
lattice sites M. We study the Hartree-Fock wave function built up from localized wave functions 
w{r) of single atoms, with nearest neighboring overlap. The zero-momentum particle number is 
expressed in terms of permanents of matrices. In one dimension, it is analytically calculated to be 
aN{M — N + 1)/M, with 01 — \ ^ w{v)dQ.'^ /[{I + 2a)l], where a is the nearest-neighboring overlap, 
I is the lattice constant, a is of the order of 1. The result indicates that the condensate fraction 
is proportional to and of the same order of magnitude as that of the vacancy concentration, hence 
there is off-diagonal long-range order or Bose-Einstein condensation of atoms when the number of 
vacancies M — N is a. finite fraction of the number of the lattice sites M. 

PACS numbers: 67.80.bd, 67.80.-s, 05.30.Jp 



I. INTRODUCTION 

Supersolidity refers to the superfluid-like behavior of 
a solid, in particular, the non-classical rotational inertia 
(NCRI) or missing moment of inertia, as a consequence of 
Bose-Einstein condensation (BEC) or off-diagonal long- 
range order (ODLRO)^'^'^. A few years ago, Kim and 
Chan observed NCRI in bulk solid ''He in torsional os- 
cillators^, which was subsequently confirmed by several 
other experimental groups^sii^i^iiS. Heat capacity ex- 
hibits a peak near the onset of NCRI". Superfluid-like 
mass flow was seen close to the melting temperature^^, 
and on melting curved, being carried by liquid regions 
at the interfaceii. Recently, it has also been observed 
off the melting curve by injecting atoms from superflu- 
idsi^. Increase of shear modulus was observed at low 
temperatures^^, but similar phenomenon in solid ^He is 
not accompanied by NCRI, indicating that elastic stiff- 
ening alone cannot produce NCRI^^. 

BEC of zero-point vacancies in the ground state of 
solid Helium is the basis in some proposals of supersolid- 
ity mechanis m'^i'^i^° . Path integral Monte Carlo studies 
indeed found that solid *He is commensurate without 
BEC^^. It was argued that zero-point vacancies or in- 
terstitials are necessary for supersolidity^^. Analytical 
calculations based on insulator-like trial wave functions 
also showed that a commensurate solid cannot be a su- 
persolid^i^^iSiiS^. But on the other hand, zero point va- 
cancies are found in variational studies using Jastrow or 
Shadow wave function^, vacancy induced BEC was also 
found by using shadow wave functions^!. A recent diffu- 
sion Monte Carlo study of commensurate solid ^He using 
another trial wave function found a condensate fraction 
~ 10"'' and a superfluid fraction < 10"^-^. 

Vacancy-based mechanism was disfavored by some re- 
searchers for the reasons that ''He is believed to be com- 
mensurate while the vacancies tend to be phase sepa- 
rated because of attractio n^^i'^°'^' . Nevertheless, the in- 
teraction between vacancies may be more complicated^. 
On the other hand, disorders such as dislocations, grain 



boundary and glassiness indeed appear to be impor- 
tan t'°'''^i'^''^'^i'^^'^^ . Grain boundaries does not seem to 
be the fundamental origin of NCRI, which has also been 
observed in large crystals'^-. In considering disorder or 
glassiness, there are theories combining this aspect with 
superfluidit y^^i'^^'^^i''"''" , as well as theories without re- 
sorting to superfluidit y''^''''^ . 

With all these results, the issue whether the ground 
state of solid ^Hc is commensurate or incommensurate 
and the mechanism of NCRI in solid ^He are still open 
questions^. It is possible that intrinsic zero point va- 
cancy is the fundamental origin of supersolidity, while 
assisted by the extrinsic disorders. This possibility is con- 
sistent with the finding in simulations that the gap for 
vacancy creation can be closed under a moderate stress^. 
Most recently, Anderson put forward a Gross-Pitaevskii 
theory of dilute gas of vacancies to account for the su- 
persolidity, arguing that every pure Bose solid's ground 
state is a supersolid based on vacancies^. 

As a theoretical approach shedding light on supersolid 
mechanism, it is interesting to consider phenomenologi- 
cal trial wave functions of a quantum solid, and exam- 
ine whether they give rise to ODLRO and supersolidity. 
One of the trial wave functions is metal-like, which is a 
product of copies of the same extended single atom wave 
function, each being a superposition of localized wave 
functions at all lattice sites. This is a BEC state, even in 
the case of a perfect crystal. First studied in 1970s, this 
wave function suffers the shortcoming that the probabil- 
ity amplitude of a configuration with one particle on each 
site tends to vanish when N — > oc^*^. Recently, it was re- 
considered with multiplication of Jastrow factors, which 
suppress multiple occupancy in a same site^. However, 
it still has the shortcoming that the equality between the 
lattice site and the number of atoms is a coincidence^'. 
In the trial wave function used in the recent diffusion 
Monte Carlo study which found BEC^^, the single parti- 
cle part is replaced as a product of wave functions on all 
lattice sites, each being superposition of wave functions 
of all possible single occupations of this site. 
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Another trial wave function is insulator-like, with 
the single-particle part being a symmetrized product of 
the localized single-atom wave functions. There is no 
ODLRO in such a wave function, even though there 
is wave function overlap between nearest neighboring 
atomaSi. The nonexistence of ODLRO was further 
proved in the cases that two particles cannot come too 
closed, and that the sum of overlap integrals of a sin- 
gle atom wave function with its neighboring ones is less 
than unitj*2^. Recently, the nonexistence of ODLRO or 
NCRI was generally shown for the case that the overlap 
between the neighboring atoms decays exponentially or 
faster, with the decay constant much smaller than the 
system size, with or without Jastrow factors^. 

On the basis of the insulator-like wave function, Imry 
and Schwartz introduced vacancies in the case that there 
is no overlap between single atom local wave functions^!. 
They found the zero-momentum particle number to be 



No^ N{M - N +1) 



\I^(.r)dn\^' 



(1) 



where M is the number of total lattice sites, N is the 
number of atoms, w(r) is the single atom wave function, 
f2 is the volume. M « L"^, sa Ml^, where L is the 
the number of atoms on each side of the lattice, d is the 
dimension, I is the lattice constant. 

But nearest neighboring overlap is crucial in a quantum 
solid. Moreover, there is some inconsistency in discussing 
BEG under the assumption that there is no overlap be- 
tween neighboring atomic wave functions. If the overlap 
is zero, | / w{r)dil\'^ also becomes zero. Then ([Ij be- 
comes not useful, as A'o ~ 0. To see this clearly, one 
can fiducially assume the single-atom wave function to 
be Gaussian, as indeed used in variational calculations of 
soHd ^He^iia, i.e. 



w{r) 
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where d is the dimension of the lattice. Then 



(2) 



(3) 



The nearest neighboring overlap is 

a = y"w(r)w(r-l)dl7 = exp[-i(^)'']. (4) 

Therefore the overlap a ^ means f ^ 0, or ^ ^ L But 
then |/w(r)drj|2 ^ 0, or \ J w{v)dn\^ / 1'^ < 1. There- 
fore, it is indispensable to consider nearest neighboring 
overlap. 

In this article, we consider the Hartree-Fock wave func- 
tion of a quantum solid with vacancies in presence of 
nearest neighboring overlap of single atom wave func- 
tions. We obtain an analytical expression for the zero- 
momentum particle number TVq, in terms of permanents 
of matrices. This expression formally reduces to Eq. ([T]) 



if the nearest neighboring overlap integral a is set to be 
0. We have made the analytical calculation of iVp in 
one dimension. Our result on A^'o indicates that there is 
ODLRO when the number of vacancies is a finite frac- 
tion of the number of lattice sites, in presence of nearest 
neighboring overlap between single-atom wave functions. 

The rest of this article is organized as follows. In 
Sec. II, we consider the Hartree-Fock wave function for 
a quantum solid with vacancies, constructed in terms of 
localized single-atom wave functions. We obtain the an- 
alytical expression for the zero-momentum particle num- 
ber, which is expressed in terms of the permanents of ma- 
trices. In Sec. Ill, we make a calculation in the case that 
the overlap integral between neighboring atoms is zero, 
reproducing the formula obtained by Imry and Schwartz. 
In Sec. Ill, we make the calculation for the case that the 
overlap integral is nonzero. Part of the mathematical 
derivation is presented in the Appendix. The summary 
and discussions are made in Sec. IV. 



II. TRIAL WAVE FUNCTION AND THE 
EXPRESSION FOR ZERO-MOMENTUM 
PARTICLE NUMBER IN TERMS OF 
PERMANENTS OF MATRICES 

We consider the following Hartree-Fock wave function 
of a bosonic solid with vacancies, 

N 

/ p, 1=1 

where w is the localized single-atom wave function, which 
is real and nonnegative, / represents a selection of N sites 
{R/j • • • R/jv} from the total M sites, Pi represents the 
A^! permutations of these selected N sites, the summation 
over / represents M\/N\{M — N) \ different choices of the 
N sites. The normalization constant A is obtained as 

= E E E E n QiPii^u) Pr (R/.')], (6) 



where 



Q(R - R') = j w{r - R)ui(r - R')dr2. 



In our consideration. 



Q(R-R') = 



1 : 


if 


R 


R' 


= 0, 


a : 


if 


|R 


-R' 




: 


if 


|R 


-R' 


>l. 



A ^ can be rewritten as 



(7) 



(8) 



(9) 



where 7'[A(/,/')] is the permanent of an A^ x subma- 
trix A(/, /') of the M x M matrix Q, whose elements 
are 



(3(R'i — Rj") 
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where Ri and Rj run over all the lattice sites. The sub- 
matrix A(/,/') is formed by choosing, from Q, N rows 
according to the set / and N columns according to the 

set r. 

The permanent of an iV x matrix A is defined as 
ViA)^ ^ e»--^"AH,-..A^.„, (10) 

where e'l = 1 when every two indices are different 
from each other, otherwise e'l ■ = 0. In other words, 
a permanent is like a determinant, except that all the 



terms in the expansion are positive, rather than with a 
sign alternation. 

In this article, the calculation is limited to a one- 
dimensional lattice. The lattice sites are numbered from 
left to right as 1, 2, • • • , AI. Q is trigonal with Qu — 1, 
Qi,i+i = Qi+i.i = a, while the other elements are 0. 
Therefore Q is 

Q = Am; 

here we introduce a square matrix A„, n being a positive 
integer, written schematically as 



/ 1 a 

a 1 a 

a 1 a 

a 1 a 



a 









V 



\ 













1 

a 






a 

1 a 
a 1 a 

a 1 / 



(11) 



where the subscript n x n indicates that it is an n x n matrix. 

For the many-body trial wave function ([5]), the one-particle reduced density matrix is 

p(r,r') ^ N J ^{r,r2, ■ ■ ■ ,rN)-9{r' ,r2, ■ ■ ■ ,rN)dfl2 ■ ■ ■ dflN 

/ /' Pi Pj, 
I I' 



(12) 



where the summation over i and j run over the rows 
and columns of submatrix A(/, /'), R; = P/(R/j), Rj = 
P/'(R/'), for which Q(Ri — R^) is just the (i,j)-th ele- 
ment Aij of the submatrix A(/, /'), Wij is the minor of 
Aij. Here the minor Wij of Aij is defined as the per- 
manent of the submatrix of A obtained by removing the 
j-th row and the j-th column. 

The number of particles at the zero momentum state 
is thus 



Nn = 



^ J p{r,r')dndn' 
XN{Q)\Jw{r)dn\'' 



YNiQ) 



n 



where 



(13) 
(14) 



y^(Q)^^^P[A(/,/')]. 



(15) 



Xn^Q) is the summation of the permanants of the minors 
of all the elements of all the N x N submatrices A's of 
Q- Yn{Q) is the summation of the permanents of all the 
N X N submatrices A's of Q. 



III. THE CASE WITHOUT NEAREST 
NEIGHBORING OVERLAP 

First let us reconsider the case without overlap between 
neighboring single-atom wave functions, i.e. a — 0, and 
show that Eq. formally reduces to Eq. ^ . 

In this case, 



Q 
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where Im represents the M x M unit matrix. Thus in 
obtaining a submatrix A(/,/'), once N rows are chosen, 
there is only one choice of N columns to give rise to a 
nonvanishing permanent. Namely, the ordering numbers, 
in the parent matrix , of the chosen columns must be 
equal to those of the chosen rows, i.e. A(/, /') must be 
a unit matrix in order to have nonvanishing permanent. 
Consequently, 

M! 

= Nl{M^N)V 

which is just the number of ways of choosing N rows. 
Note that the order of the chosen rows and the order 
of the chosen columns both remain the same as in the 
parent matrix. 

In order to calculate Xn{Xm), we need to find out 
all nonzero minors for all submatrices A's of Tm- Note 



that 'P(A) = does not mean A has no nonzero mi- 
nors. Given that the parent matrix is a unit matrix 
in order that A has one or more nonzero minors, 
A must be either a unit matrix Xjv or diagonal with only 
one "0" diagonal element. In the former case, there are 
Ml/Nl{M — N)l ways of making up the unit submatrix 
A = Xjv, which has N nonzero minors, each equal to 1. 
In the latter case, one first choose — 1 rows and — 1 
columns, with the same ordering numbers in the parent 
matrix Xjv , to make up TV — 1 diagonal elements "1" . The 
number of ways of doing this is M!/(iV-l)!(Af-iV+l)!. 
To choose the remaining one row and one column, their 
ordering numbers in the parent matrix X^r must be differ- 
ent, such that the remaining diagonal element in A is "0" . 
The number of ways of doing this is {M — N + 1){M — N). 
Each A so obtained only has one nonzero minor, which 
is equal to 1. Therefore, 



M 



Ml. AT 

Ar!(M-iV)!-^(^^ 



M! 

{N-iy.{M-N+iy. 

N + 1). 



{M - N){M - N + 1) 



(17) 



r 



Substituting ^ and ^ into Eq. ^ indeed recov- 
ers Eq. H]). 



IV. THE CASE WITH NEAREST 
NEIGHBORING OVERLAP 

With nearest neighboring overlap, the zero-momentum 
particle number is 



No = 



XN{AM)\Jw{r)dn\' 
Yn{A 



Ml 



(18) 



where the matrix Am is as defined in Eq. (fTTj) . 

For an arbitrary matrix Smxn, we introduce Yk{Smxn) 
and XkiSmxn), with k < min(m, n). Yk{S 

mxn) is the 

sum of the permanents of all the k x k submatrices of 
Sm X n ■ Xk{Smxn) is the sum of the permanents of aU 
the minors of all the k x k submatrices of Smxn- 
First, we note the existence of the relation 



XkiA,) = (n-fc + l)2n._i(A„), 



(19) 



for the following reason. Every minor of a fc x fc submatrix 
of An is in fact a (fc — 1) x (fc — 1) submatrix of An, while 



J 



a (fc — 1) X (fc— 1) submatrix is a minor of many different 
fc X fc submatrices. For a given (fc — 1) x (fc — 1) submatrix 
of An, one can add an additional row and an additional 
column of the An, making up a fc x fc submatrix of An, 
of which the concerned (fc — 1) x (fc — 1) submatrix of An 
is a minor. There are {n—k + 1)^ ways to do this. Hence 
a (fc — 1) X (fc — 1) submatrix is a minor of (n — fc + 1)^ 
different k x k submatrices of An, thus we obtain the 
relation (fT9l). 



Therefore 



Xk{An) _ (71-fc+l)2yfe„i(A„) 



Yk{An 



Yk{An) 



(20) 



In Appendix A, we obtain that for n > 2, 



Y^iAn) =(l + 2a)' 



, n{n — 1) 



(5a^ -I- Aa)n + 7a^ + 4a 



„ ,^n(n — 1) 



1. 



2! Il + °<^'1- 



where 0{l/n) represents a term of the order of 1/n. 
In Appendix B, we obtain that for any 3 < fc < n. 



(21) 



k—2 n—s—1 



n-2 



Yk{An) = J2 Yk-iiAi) + 2Y,o.' Yk-siAi)+a^ ^ ^-2(^0 

l = k-l s=l l=k-s 

+ (1 + 2a)2a''-\n ~k){n-k+ 1). 



l = k-2 



(22) 



5 



From this relation, we know that Yfc(A„) > increase of s. Since a < 1, RHS of ([22)1 converges with 

Yfe_i(A„_i), yfc_i(A„_i) being merely one term in the respect to a. Also note that the last term is of the power 

first summation in RHS of (|22p . Consequently, in the of a'^"^. 

summation over s, '2YM=k-s ^fc-s(^i); which also de- 
pends on a, by which a* is multiplied, decreases with the Therefore, Yk(An) can be written as 



n-2 



n(A„) = [J2 Yk-i{Ai) + 2a J2 Yk-iiAi)][l + 0<ia)] 

l=k-l l=k-l 

= [(l + 2a) ^ Yk-i{Ai) - 2aYk-i{A„^i)][l + 0<{a)], 



l=k-l 



r 



(23) 



where 0<{a) denotes a term at most of the order of a. with k — 1 < I < n, i.e., 
In the following, we show by induction that 



y,(A„) = (1 + 2a) \ [1 + 0<{a)]. (24) ^ + (fc - 1)!(Z fc + 1)! 



k\{n-ky.' 

Suppose that the similar identity is valid for Yk-i{Ai), Then 
I 



f[l + 0<(a)]. 

(25) 



(l + 2a) J2 Yk-i{Ai) 



k 



l=k-l 



(l + 2a) 



E 



[l + 0<(a)]. 



(26) 



Using the identity^ 
p 

we obtain 



Hence (1^ becomes 



1 (p + g + 1)! 
9 + 2 (p-1)! ^ 



E 



1 n\ 



l=k 



^^{l-k + iy. k{n-k)V 



(27) 



(28) 



(l + 2a) ^ r,_,(A0 = (l + 2a)'=^^^-!^[l + O<(a 



(29) 

On the other hand, according to the assumption ((25|) 
and the above result ([29]) . we have 



2arfc_i(A„_i) = 2a(l + 2a) 



fe-i 



(n-1) 



2a fc 



(fc- l)!(n-fc)! 

ri-l 



[l + 0<(a)] 



l + 2an 



[(l + 2a) ^ y,_i(^0][l + O<(a)] 



l=k-l 



0{a)- Yk^Ai), 



(30) 



where k/n < 1. 



Substituting (HH) and ([201) into Eq. ([121) yields Eq. 
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hence the proof completes. 

Then we substitute the proved identity into 
Eq. ([20]), obtaining 



Xk{A^) k{n-k+l) 



1 + 2a 



(31) 



if a < 1. 

Therefore, according to p^ . the zero- momentum par- 
ticle number is 

iVo « N{M-N + l)\Sw{v)dn\-' ^ ^^2) 



1 + 2a 



This identity formally reduces to that obtained by Imry 
and Schwartz long ago when we set a = 0. But 
|/w(r)dJ7|2 ^ only if a 7^ 0. 

On the other hand, in case 71 3> A;, no matter whether 
a <C 1, we can also obtain the identity pip and thus the 
result ([32]) . This condition does not correspond to the 
physical situation concerning solid *He, as that would 
mean most of the lattice sites are empty. For complete- 
ness, we give the mathematical proof for this case in Ap- 
pendix C. 



0.29 for solid 



,49 



usmg 



/372e/Z 



under the Gaussion wave function ([2]). For (i = 3, l/(l-f 
2a) should be replaced by another function /(a) of a, 
which should still be of the order of 1. Anyway, a = 
{2y/T:£_/lY f{a) must be of the order of 1. 

Therefore there is BEG of atoms, i.e. A^o is a finite 
fraction of A, when the number of vacancies M — A^ is 
a finite fraction of the number of lattice sites M . This 
condition also implies that the number of atoms A is a 
finite fraction of M . 

Interestingly, the condensate fraction Nq/N is pro- 
portional to and of the order of vacancy concentration 

{M - N) /M, 



No 
A 



M -N 



(36) 



Gurrently, the experimental upper bound of vacancy con- 
centration is about 0.4%2^. Hence a Hartree-Fock wave 
function for a solid with zero point vacancy implies that 
the condensate fraction is about 0.004q;, which is very 
reasonable. 

Moreover, for such low vacancy concentration, one has 



V. SUMMARY AND DISCUSSIONS 

To summarize, we have studied the Hartree-Fock wave 
function of a lattice of atoms with N < M, where M 
and A^ are the numbers of lattice sites and atoms, respec- 
tively. The Hartree-Fock wave function is constructed in 
terms of localized wave functions of single atoms, with 
nearest neighboring overlap. 

In one dimension, under this wave function, we have 
obtained the zero-momentum particle number as given 
in (1321). which can be rewritten as 



A^( 







A^o = a 



where 



A^(Af-A^+l) 
M ' 



(l + 2a)/ 



(33) 



(34) 



is a finite fraction of the order of 1. 

To be specific, let us again use the Gaussian wave func- 
tion for w{r), as given in ([2|. Then in one dimension, 



20FS 



[l + 2cxp(-i(|)2]' 



(35) 



which is of the order of 1 when ^/l is a finite fraction 
around 0.36, which is obtained from the Lindemann ratio 



M - N 



(37) 



i.e. a equals the number of condensed atoms per vacancy. 
This is well consistent with the result of variational sim- 
ulation based on Shadow wave function, which gives 0.23 
condensed atoms per vacancy at 54 bar^^. 

The Hartree-Fock wave function could be the ground 
state of a mean field theory. Although it is not multiplied 
by the Jastrow factor, the double occupancy is excluded 
by construction. Our calculation is done for one dimen- 
sion. In three dimensions, there should not be qualitative 
difference in order of magnitude from the result for one 
dimension. Hence our result is qualitatively informative 
for solid ■^He, suggesting that its supersolidity based on 
BEG of atoms induced by zero point vacancy is possible. 
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APPENDIX A: YiiA^) 



Here we calculate Y2{An) {n > 2), i.e. the sum of the permanents of all the 2x2 submatrices of An- 
In addition to the definition of j4„ as given in (fTT|) . we shall also use matrices 



/a a 

1 a 

a 1 a 

a 1 









Vo 



C'(„-2)x(n-l) 



/ a 1 a 



a 



a 1 

a 1 

a 

■■• 





\ 



\ 











a 

1 a 

a 1 a 

a 1 / 



(n-l)x(n-l) 



\ 





a 





■. 1 a 

a 1 a 
a 1 / 







(n-2)x(n-l) 



The permanent of an n x n matrix S is equal to 



where Wij is the minor of the Sij 



Using this property, we can expand l2(^n) as 



/ 1 a 

a 1 a 

a 1 a 

a 1 

a 



\ 











a 







V 



■•• ■•• 

a 1 a 
a 1 a 
a 1 a 
a 1 / 

= Yi(A„_i) + aYi{B(n-l)x{n-l)) + aYi{C(^n-2)x(n-l)) + i1j(A„_i) 

= Y2{Ar,-i) + YiiAr,-i) + 2ayi(A„_2) + 3a2, 

where we have used yi(S(„_i)x(„-i)) = 2a + Yi(A„_2), ^i(C(,i-2)x(n-i)) = a + Yi(A„_2), which are straig 
In this way, we obtain the following set of identities 

F2(A„) - Y2{Ar,-i) = Yi(A„_i) + 2ari(A(„_2)) + Sa^, 
y2(A„_i) - y2(A„_2) = Yi{An-2) + 2ari(A(„_3)) + 3a2, 



Y2{A3) ~ Y2{A2) = Yi{A2) + 2aYi{Ai) + 3a\ 
Adding these identities together gives rise to 

n-2 

>2(A„) - Y2{A2) = ri(A„_i) + (1 + 2a) J2 ^i(^j) + 2an(Ai) + 5a^{n - 2). 



8 



Clearly 12(^2) = 1 + a^, Yi{Aj) = j + 2{j — l)a. Hence it can be obtained that 



Y2{An) = (1 + 2af - (5a2 + Aa)n + Ic? + 4a, 



(A7) 



which is also satisfied when n = 2, as 1^2(^2) = 1 + 



APPENDIX B: Yfc(A„) 



We now calculate Yfe(A„) for 3 < fc < n, in a way similar to the calculation of y2(^n) above. 
Similar to (jA4[) . we obtain 



/laOOO 000\ 

alaOO 0000 

OalaO 0000 

OOala 0000 



a 







0000a •• -.000 

00000 a laOO 

00000 alaO 

00000 Oala 



VOOOOO OOal/ 

Yk-i{An-i) + aYfc_i(i?(„_i)x(„-i)) + ayfe_i(C(„_2)x(ri-i)) + ^felAi-i) 



(Bl) 



Expansion of Yfe_i(B(„_i)x(„_i)) gives 



Yfe_i(B(„_i)x(„-l)) aYfc_2(^fc-2) + ^fe-l(C(Ti-l)x(n-2))- 

For any matrix S, Yk-i{S'^) — Yk^i{S). Hence we have 

Yk{A,) = rfe(A„_i) + Yk-i{An-i) + a^Yk-2{Ar,-2) + 2arfe_i(C(„_2)x(«-i))- 



(B2) 



(B3) 



Iterative expansion of yfe-i(C(„_2)x(n-i)) yields 



5^fc-l(C(„_2)x(ri-l)) — Yk-l{An-2) + aYk-2{C(n-3)x{n-2)) 



fe-2 



a'' ^Yk^s{A„-s-i) + a'' ^Yi{C(^n-k)x{n-k+i)), 



(B4) 



where yi(C(„_fe)x(„-fe+i)) = (n - k) + a + 2{n - k - l)a. 
Therefore 



rfe(A„)-rfc(A„_i) 



yfc_l(A„_i) + a^Yk-2{An-2) + 2ayfc_i(C(„_2)x(n-l)) 

yfe_i(A„_i) + a^Yk-2{An-2) 

fe-2 

+2^a''Yk-s{An-s-i) + 2a''-^[{n -k)+a + 2(n -k~ l)a 



(B5) 



Hence 



fe-2 



+2 ^ a"yfe_, (A„_5_i) + 2a'''"i [a + (n - fc) + 2(n - fc - l)a] 

fe-2 

+2 ^ a^yfe_, (A„_,_2) + 2a^-i [a + (n - 1 - fc) + 2(n - A: - 2)a] 



s=l 



>fc(Afc+2) - >fe(^fe+i) = Yk-i{Ak+i) + a^rfc_2(Afc) 

fc-2 

+2Y,ci'Yk^s{Ak+i-s) + 2a'=-i[a + 2 + 2a]} 
Yk{Ak+i) - Yk{Au) = n/r(Afe) + a'Yk-2{Ak-i) 

fc-2 



+2^ a^yfc_,(Afc_,) + 2a^-i[a + 1]} 



Adding these identities together leads to 

n— 1 n — 2 /c — 2 n — s — 1 



(B6) 



(B7) 



ifc(A„)-yfc(Afc) + E n-2(A,) + 2^a^ ^k-s{Ai) 

l—k l—k—1 s—1 l—k—s 

n—k— 1 

+2a'=-i ^ [a+(n-fc-s) + 2(n-/c-l-s)a] 

s=0 

Since Yml^n) exists only when m < n, we have 

>fc(A-) = n-i(A-i) + a'i"fc-2(Afc-2). (B8) 

Therefore 

n— 1 fc— 2 n— s— 1 n— 2 

YkiA^) - ^ n._i(A0 + 2^a^ ^ Yk-s{Ai)+a^ ^ rfc-2(Az) 

i=fc-l s=l l=k-s l = k-2 

+ {l + 2a)2a^~^{n-k){n-k + l). (B9) 



APPENDIX C: CALCULATION OF No IN THE CASE OF ti » fc 

Here we show that the identity (|3ip and thus the result ([5^ are also valid if n ^ fc, no matter whether a <C 1 or 
not. Mathematically, n S> fc means n — s- cx) while fc remains finite. 
In the following, we show by induction that 

yfc(A„) = (1 + 2a)'=^j^^[l + 0{a)0{\)]. (CI) 

Suppose that the similar identity is valid for Yk-s{Ai), with s > 1 and k — s < I < n, i.e., 

y._.(AO ^ (1 + 2a)-^^^— + OHO(i)]. (C2) 

In the second term in the exact identity (|B9[) for yj,(A„), 2a* is multiplied by Y^^i^l^I^ Yk_g{Ai)^ which can be 
evaluated by using the assumption (|C2|) to be 

l=k-s ^ ' l=k-s ^ ' 
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(l + 2a)'=-" {n~-s)\ 
{k-s + iy. (n-fc-1)! 

f V y (A ^l (fc-2)...(fc-,s + l) 1 



where we have used the identity ([27] 

nil I T" 1 -rx 1 1 ^ /H r \ 



In the third term in (|B9[) . 2a^ is multiphed by jyi=k-2 ^k-2{A-i), which can be similarly evaluated to be 



/=fc-2 ^ ^ l=k-2 ^ ' 

(l + 2a)*--2 (n- 1)! 
(fc-1)! {n-k)\ 

■n-2 , 

where we have also used the identity (I27p . 

Besides, the last term in is 0(a''"i)0(n2) < Yfc_i(A„_i) = ©(n'') if fc < n. 
Therefore 

ri-l n-2 

yfc(A„) = yfc_i(A0 + 2a i^fc-i(^0][l + O(a)O(-)], (C5) 

where 0(l/n) ^ 1, 0(a) is of the order of a, which we do not need to specify. 
Using the assumption (|C2p for s = 1, we obtain 



where 



71 — 1 n — 1 T| 

(l + 2a) E r._,(A.) ^ (l + 2a)'= ^ (fc , ^ i), + ^(7)] 



Using the identity (P7)l . we obtain 



Hence (IC6I) becomes 



(l + 2a) I; n-.(/l,)= "''J°''(„';'^, [l + 0(a)0(l)]. (Cll) 



11 



On the other hand, according to (IC2p and (|C11|) . 



for the reason that k <^n. 
Therefore, 



2a k, 



71 — 1 ^ 

[(l + 2a) Yk-i{Ai)]0{a)0{-), (C12) 



n 

=k-l 



Yk-i{Ai)+2a Y Yk-i{Ai) 

l=k-l l=k~l 
n-1 

(l + 2a) Y Yk-iiAi) - 2aYu-i{Ar,-i) 

n— 1 

= (l + 2a) V Yu-i{Ai)[l + 0{a)0{-)]. (C13) 



By using (jClip and (|C13[) . (|C5[) leads to the identity (jCl|) . which is thus proved. Substituting this proved identity 
into Eq. (gOl), we obtain Eq. and thus also ((5^ . 

I 
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